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Effective star formation rates in tabular form are computed which yield a prescription for 
the star formation activity in model galaxies as a function of ambient density, metallicity, and 
stellar feedback. The effects of supernova explosions on the thermal balance of the Interstellar 
Medium (ISM) and the presence of a multi-phase ISM are explicitly included. The resulting grid 
, of models can be implemented easily in N-body codes for the computation of star formation 

processes in merging galaxies and cosmological simulations. 



subject headings: galaxies: star formation - galaxies: structure - ISM - molecular processes 



Introduction 



> 

m 

<N 
O 
OO 

g> 

i-^ ' The past decade has seen tremendous advances in the field of numerical cosmology (Navarro, Frenk & 

^P^i White 1996 and references therein). N-body simulations and hydro codes have benefitted from the large 

increase in CPU and memory capabilities of the current generation of supercomputers, and have reached 
' a high enough sophistication to describe the development of large scale structure in the universe and the 

formation of galaxies with redshift. An important issue which ultimately needs to be addressed, in order 
to compare theory with observations of high redshift galaxies, is the formation of stars in proto-galactic 
structures (Norman & Spaans 1997, hereafter NS97, and references therein). Many parameterizations of 
the star formation rate in the ambient Interstellar Medium (ISM) exist in the literature and mostly derive 
from a Schmidt law applied to a sufficiently large body of gas (Theis et al. 1992; Kauffmann & White 1993; 
Spaans & Norman 1997, hereafter SN97). 

Star formation is a local phenomenon which must find its explanation in the stability and fragmentation 
of dense molecular clouds. Studies in our own Galaxy have focussed on the structure of dense proto-stellar 
cores, along with the chemical and thermal balance of star-forming regions (Helmich 1996, and reference 
therein). These studies lend indirect support to a Schmidt law, but emphasize the need to include explicitly 
the structure of the multi-phase ISM to model accurately the most important heating and cooling processes. 
A large unknown in these investigations is the relative importance of feedback. Supernova explosions and 
stellar radiation associated with the process of star formation influence the global physical structure of the 
interstellar gas which supports this process. 

Ideally, in the cosmological context one would like to solve for the properties of the thermal and 
chemical balance of the ISM and the star formation rate simultaneously with the solution of the gravitational 
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N-body and hydrodynamical problem. Such an approach is beyond what is currently feasible and one is 
forced to decouple the ISM and star formation problem from the dynamical one. Along these lines, the 
purpose of this work is to present a set of numerical simulations which span a large range in density and 
metallicity, include a quantitative prescription for feedback, and yield the resulting star formation rate in 
the ambient medium in tabular form. This work extends the methods developed in a series of papers where 
the structure of the ISM and the importance of feedback effects have been investigated for proto-galactic 
disks (NS97), dwarf galaxies (SN97), massive proto-spheroidals in the Hubble Deep Field (Spaans & Carollo 
1997), and elliptical galaxies (Carollo & Spaans 1998). 

The strength of the procedure used in this work to produce the effective star formation rates is its 
detailed treatment of molecular line cooling, chemistry, and thermal balance (section 2) for a wide range of 
physical conditions. This allows for a detailed treatment of the multi-phase structure of the ambient ISM, 
and the feedback of stellar photons and supernova ejecta on its physical balance. The star formation rate is 
computed in our approach through a Schmidt law applied to the molecular phase of the ISM. The use of 
this empirical law presumes that star formation as observed in our local environment is representative of 
star formation in general. This seems a reasonable approach for spiral galaxies (Kcnnicutt 1989). 

In the following sections, the input physics are discussed. The processes involved are complex, and 
various physical phenomena influence the final outcome of the computation. Therefore, one should keep 
the following major calculational steps in mind. 1) The thermal and chemical balance of the medium is 
determined locally and used to compute the amount of molecular gas. 2) This allows a calculation of the 
formation rate of stars according to a Schmidt law and an initial mass function (IMF). 3) The produced 
radiation is propagated across the computational grid in a radiative transfer calculation, and used to 
determine heating, ionization and dissociation rates. 4) The supernova blast waves, produced as the end 
products of stellar evolution, are propagated across the grid as well and yield the local input of kinetic (gas 
bulk motions) and thermal (gas heating) energy, as well as metals. 5) Steps 3 and 4 provide input for step 
1. This entire procedure is followed for many different ambient densities, to derive the enrichment and 
physical state of the gas together with the corresponding star formation rate. 



2. The Model: Stars Embedded in a Multi-Phase ISM 

The models include three stellar and three gaseous components. The gas phases include the cold 
molecular clouds, the warm neutral/ionized medium, and the hot tenuous interiors of supernova bubbles. 
The phases are assumed to be in pressure equilibrium and their chemical and thermal balance is computed 
explicitly. The theoretical background for the evolution of the multi-phase ISM in primordial galactic 
structures is described in NS97 and SN97, where the employed numerical methods and a discussion of 
the spatial and temporal resolution of the adopted grid can also be found. Although no dynamics are 
included explicitly, there is mass exchange between the cold molecular and hot tenuous phase due to cloud 
evaporation. The density dependence for mass evaporation is given by oc E 6 / 5 n^ 4 ' 5 , with E the supernova 
energy and rih the density of the hot phase (McKee & Ostriker 1977). In the following we adopt for a grid 
cell a typical size of 20 pc, and time steps of 30-60 Myr. Even though we determine stationary solutions, it 
is necessary to integrate in time to determine the response of the ambient gas to the star formation process, 
i.e. stellar evolution, propagation of supernova shock waves etc. 

The stellar components are divided according to their final evolutionary stages into massive stars 
of more than 11 M©, which explode as Type II supernovae, and low-mass stars, which are assumed to 
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loose their material in a planetary nebula phase. The low-mass stars will become 0.6 M white dwarfs 
(Weidemann & Koester 1983). Stars with masses below 0.6 M Q do not evolve during the lifetime of a 
galaxy. The third class of stars comprises the stellar remnants in the form of white dwarfs, neutron stars, 
and black holes. A fraction of the white dwarfs gives rise to Type la supernovae in the merging CO-dwarf 
picture (see below). Below we describe in detail the set of equations which are solved for the star formation 
(SF) and feedback processes: chemical and thermal balance, the occurrence and propagation of supernovae, 
metal yields and enrichment, radiative transfer, and stellar life cycles. 



2.1. Star Formation 



Following Larson (1991) the star formation rate is calculated with a Schmidt (1959) law applied to 
some volume in a galaxy with a gas mass M 

/ ^cm \ h / 1 \ 

= osf(-o-), (1) 
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where the label "cm" indicates the cold molecular phase, n cm is a number density, n° m = 40 cm 3 and 
b = 1 — 2. The coefficient asF is normalized to the SF rate as observed in molecular gas in the Solar vicinity 

M 

«sf = -g- M yr \ (2) 

with M normalized to the total body of gas M g in the Milky Way (see below). For our own Galaxy as a 
whole this corresponds to an instantaneous rate of ~ 3 M Q yr _1 (Wyse 1997, private communication). The 
generic value of b is 1.3, but b = 2 is adopted in those low metallicity regions where the ambipolar diffusion 
time is shorter than the free-fall time of a molecular cloud (see NS97 for a detailed discussion). 

The density of the cold molecular phase follows from application of the Field, Goldsmith, and Habing 
(1969) formalism to the thermal balance of a multi-phase ISM. The angular momentum of star-forming 
clouds is not included beyond the accuracy of the local Schmidt law applied to the cold molecular 
component. As such, we cannot address the (important) question of cloud fragmentation, and the IMF is 
an input parameter. We have assumed a time- independent Salpeter IMF <j>{m) oc to~ 2,35 . Stellar masses 
between 0.1 M and 40 M are considered. It follows that approximately 10% of the newly formed stellar 
mass is incorporated into massive stars with a lifetime of no more 2 x 10 7 years. 



2.2. Chemical Evolution 



The evolutionary end stages of stars produce metals which provide most of the cooling (atomic and 
molecular lines) and dust opacity in a model galaxy. It is thus vital to know the production and spatial 
distribution of metallic species. 



2.2.1. Elemental Abundances in the Gas Phase 
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The chemical evolution of elemental species has been discussed extensively in the literature, and 
the equations as described in Timmes, Woosley & Weaver (1995) and Tantalo et al. (1996) are adopted. 
The main difference with our implementation is the explicit three-dimensional spatial dependence of the 
chemical quantities. The treatment of physical processes is explicitly three-dimensional, and the spatial 
enrichment of the model galaxy through the expansion fronts of multiple supernovae is included. While 
the supernova bubble expands with time, the rate at which metals condense out into the cold molecular 
phase is computed to determine the local enrichment. The cooling rate of the enriched gas then determines 
whether a molecular phase can be supported, and at what interstellar gas density the SF process takes 
place. This eliminates the assumption of a one zone system. Allthough there is no infall of material, and 
the model galaxy is therefore a closed box, there is mass exchange (in particular metals) between the grid 
cells through the action of the SN expansion fronts. 

For the stellar yields, three distinct components are specified: one for the massive stars 
(11M Q < M < 40M Q ) which become Type II supernovae, one for the intermediate to low mass 
stars (M < HMq) which become planetary nebulae, and one for the remnants of the intermediate to 
low mass stars that become Type la supernovae (Timmes et al. 1995). The Type II supernova yields of 
Woosley & Weaver (1995) are adopted throughout. The yields from intermediate to low mass stars were 
taken from the CNO models with 1M < M < 8M of Renzini & Voli (1981). We have adopted their case 
B table for Z — O.OCMZq and Z — Zq, and applied linear interpolation for the intermediate metallicities. 
For metallicities larger than Solar we extrapolated the yields linearly. This extrapolation does not strongly 
influence our results. By the time the ISM has been enriched to super-Solar metalliities, most of the baryons 
in a model galaxy have been incorporated into stars. The W7 model of Nomoto, Thielemann & Yokoi 
(1984) was adopted for the Type la supernovae. 

The average mass £,i(p, t) and abundance Xi(p, t) of the elemental species in a grid cell centered on the 
point p — (x, y, z) at time t are governed by 



*MgQ=-X i (p,t)S(p,t)+ r B S(p,t-t M )Y M ,i(t-tM)I(M)dM+ 

r-M B /-0.5 

C I(M)[ f( f i)S(p,t-t M2 )Y M ,i(t-tM 2 W}dM+ 

JlTlB J Mmin 



r M B 

(1 - C) / S(p, t - t M )Y M ,i(t - t M )I(M)dM+ 

J ran 



/ S(p,t-t M )Y M ,i(t-tM)I(M)dM. (3) 

J M B 

In the equation, S(p, t) is the normalized SF rate, Ym,% are the elemental yields of elements i from stars 
of mass M, I(M) is the initial mass function whose lower and upper limits are M m { n and M max , tM is 
the lifetime of a star of mass M, and £m 2 is the lifetime of the secondary in the binary, which determines 
the time required for the system to undergo a Type la event. In the models this lifetime depends on the 
initial chemical composition as tabulated in Bertelli et al. (1994). The various integrals which appear in 
Equation (3) represent the contributions of Type II and Type la supernovae in the formulation of Matteucci 
& Greggio (1986). Specifically, the second integral reflects those binary systems which have the necessary 
properties to produce Type la supernovae, tub and Mb are the lower and upper mass limit for the total 
mass of the binary system, f(p) is the distribution function of their mass ratios, and f-im'm 

denotes the 

minimum value of this ratio. The lower integration limit is equal to 3 M Q , to ensure that the accreting 
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white dwarf can reach the Chandrasekhar limit. The upper mass limit is equal to 16 M Q for a binary 
system, based on the assumption that the maximum mass of the primary that produces a carbon-oxygen 
white dwarf is 8 M©. The constant C indicates which fraction of all binary systems undergo supernova 
events. This constant is chosen equal to C = 0.007 to reproduce the solar 56 Fe abundance (Timmes ct 
al. 1995). 



2.2.2. Atomic and Molecular Species 



If one knows the local elemental abundances and the radiation field (discussed below) then it is 
possible to determine the chemical composition of the ambient gas. To derive the abundances of atomic 
and molecular species, the solution vector of a chemical network in point p is computed according to 

(T( P W P )) = = 0. (4) 

Here H denotes the chemical network, T is the ambient temperature, and the chemical abundances by 
number m are determined in steady state. That is, the evolution time of a galaxy, a free-fall time <~ 500 
Myr, is assumed to be long compared to the time to reach chemical equilibrium, <~ 50 Myr for diffuse H 2 gas 
and much shorter for all other molecular species. Only in a very low metallicity environment, Z ~ 0.005Z©, 
is the formation time scale of H2 as long as 500 Myr and intimately linked to the dynamical evolution of 
the system. 

The adopted reaction network is based on that used by van Dishoeck & Black (1986, 1989) to model 
the chemistry of translucent clouds, combined with the sulfur chemistry as described in Drdla et al. (1989). 
The network has been checked against the UMIST database (Millar et al. 1991), and only minor differences 
have been found. The chemical network includes 24 elements, and isotopes of H, C and O. The inclusion 
of many metals is important because they dominate the ionization balance and differences in ionization 
potential become important for extinctions of more than 2 mag. The network includes 1549 reactions among 
215 species. We adopt the line and continuum photo-dissociation cross sections based on the literature 
summarized by van Dishoeck (1988). 



2.3. Thermal Balance 



The heating of the ISM involves radiative sources like stars as well as mechanical sources like stellar 
winds and supernova expansion fronts. Cooling is provided by line emission from hydrogen and metallic 
ions as well as continuum processes. In order to determine the thermal state of the gas, one needs to find 
the solutions (p,T) of the fundamental energy (E) equation 

dE{T dt P,t) =T(T,p,t)-HT,p,t). (5) 

In this equation, T is the temperature, p is the gas density, T is the total heating rate and A the total 
cooling rate, both per unit of mass. In the following subsections the various terms which contribute to the 
overall balance are discussed. 
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2.3.1. Heating by Supernovae 

The hydrodynamics of interstellar gas is not treated explicitly in our work, and the results of various 
numerical studies are used to implement the contribution from mechanical processes such as supernova blast 
waves. We adopt the formulations of Bertelli et al. (1994) for Type II, and of Greggio & Renzini (1983) for 
Type la supernovae. The total energy injection due to supernovae of Type X which last for a time At is 
given by 

T x = [ e x (t-t')R x (t')dt', (6) 

J At 

where Rx(t) is the number of supernovae per unit mass and time, and ex(t) describes the thermalization 
of the supernova energy as a function of time when cooling effects are properly taken into account. For ex 
we adopt the results of Gibson (1995) and Thornton et al. (1997). From the latter authors we adopt the 
radial evolution of the supernova remnants to determine the distribution of metal-rich tenuous gas as it 
condenses out into the cold molecular phase. The rate Rx is computed from the local star formation rate 
in grid point p at time t and depends on the assumed IMF. We have also included in Equation (5) the work 
performed by supernova shocks through the acceleration of molecular clouds. From Thornton et al. (1997) 
we take the fraction of the initial supernova explosion energy which ends up as kinetic energy of the ISM. 



2.3.2. Heating by Stellar Winds 

The rate of energy injection by stellar winds has a similar form, namely: 

T w = e w (t-t')R w (t')dt' ', (7) 

J At 

where Rw{t) is the number of stars per unit mass and time expelling their envelopes during the time 
interval At. The quantity e\y has the same meaning as ex- A star which is shedding mass, deposits an 
energy of (Gibson 1994) 

e wo = 0.15M e ,(M)(Z/Z Q )°- 75 v(M) 2 , (8) 

where M e j is the amount of mass ejected by a star of mass M and v(M) is the velocity of the ejected 
material. The velocity v(M) is assumed to be equal to the maximum of the terminal velocity of the wind 
w 20 km s _1 or the galaxy velocity dispersion. The adopted process of thermalization is as described in 
Bertelli et al. (1994). The magnitude of this heating source is not that accurately known, but it docs not 
dominate the thermal balance of the ISM. 



2.3.3. Heating by Ultraviolet Photons 

The ionization heating and charge balance of the medium is explicitly included through the Saha 
equation applied to H, He, C, O, N, S, Fc, Si, and Mg. The energy density and shape of the radiation field 
follows from 

E vv = [ e vv (t-t')R vv (t')dt', (9) 

J At 
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where all the symbols have their usual meaning. In the determination of euv we use the stellar population 
models of Leithcrcr ct al. (1996) for a Salpeter IMF. 

To assess which fraction of the ultraviolet energy is available for heating, ionization and dissociation, a 
continuum radiative transfer calculation is required. In the literature (Granato et al. 1997 and references 
therein) one finds a typical number of 0.01 for the fraction of the ultraviolet light which is not absorbed by 
dust and re-radiated as infrared radiation. Because the abundances of metals and therefore of dust may 
vary strongly across the galaxy, the Monte Carlo radiative transfer method of Spaans (1996) is adopted to 
compute the strength of the average radiation field from the SF rate and the local metallicity (SN97). The 
star formation rate used in the radiative transfer computation is a local average over a grid cell, and yields 
a local source function for the ultraviolet radiation field. This way local fluctuations in the radiation field 
strength are taken into account in the chemical and thermal balance. 

For the molecular phase with metallicities larger then 1% of Solar, the dominant heating mechanism is 
the emission of photo-electrons from dust grains. The formalism of Bakes & Tielens (1994) is adopted with 
the implementation as described in Spaans et al. (1994). Part of the photo-electric heating of the molecular 
ISM is provided by large molecules like PAHs (polycyclic aromatic hydrocarbons), and their abundance is 
assumed to be equal to 10% of the gas phase carbon. 



2.3.4- Heating by Proto- Galactic Collapse 



Since the purpose of the effective SF rates is to study the gas phase processes during the formation 
of galaxies, we also include the heating which results from the infall of material on a time scale rg and 
leads to the formation of the galaxy (Binney & Tremaine 1987). This is an aspect of our calculations 
where a coupling with realistic gas dynamics from a cosmological formalism would be most desireablc. 
Decoupling the ISM from the dynamics, we adopt a simplified approach here. If Ri is the initial radius of 
the proto-galaxy, then 

■* = ^Sm> (10) 

where M is the total mass of the galaxy. If no other energy sources are present, a proto-galactic cloud at 
the virial temperature cools on a time scale 

t coo1 = 6.3 x 10 5 (iVl0kpc) 2 yr. (11) 

At the end of the collapse phase the system is assumed to have an energy given by 

E v i T = _E v i r . (l E22l) 5 (12) 

TS 

where .E v ir,o is the initial virial energy. The remaining energy at the end of the collapse phase is converted 
into heat at the rate 

r c = (13) 

Tff 

This source of heating remains active as long as the age of the galaxy is less than ts, and sets the initial 
gas temperature. The latter is very important for the early low metallicity chemical balance which pertains 
in the primordial gas. We adopt a collapse factor A w 30 (NS97), i.e. Ri <~ 300 kpc. 
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2.3.5. Relative Importance of Heating Processes and Galactic Winds 

As mentioned above, the heating of the molecular gas is dominated by photo-electric emission for 
metallicities in excess of O.OIZq. It is this rate which determines to a large extent the thermal and chemical 
balance of the multi-phase ISM, and therefore the effective star formation rate supported by the ambient 
molecular gas. Nevertheless, the other heating terms are of importance as well during certain epochs. The 
proto-galactic collapse phase provides the initial thermal structure of the ambient gas prior to the formation 
of stars. It thus fixes the chemical abundances of the primordial galaxy. During periods of massive star 
formation, the energy input from supernovae becomes comparable or even larger than the photo-electric 
heating rate in small regions of the galaxy where the supernova blast wave interacts with the ambient 
interstellar gas. 

If the supernova energy input dominates during some period of time, the temperature or kinetic energy 
of the interstellar gas can increase and may become comparable to the gravitational binding energy of the 
model galaxy. In this case a wind develops which is assumed to quench the SF process (e.g. Larson 1974; 
Carlberg 1984; Arimoto & Yoshii 1987; Matteucci & Tornambe 1987; De Young & Heckman 1991). 



2.4. Cooling by Line and Continuum Processes 

The cooling rate depends crucially on the abundance of certain atoms and molecules like C, O, H 2 
and CO, and the ability for radiation to escape from the medium, i.e. radiative transfer. For temperatures 
above 10 4 K, the metallicity dependent tabulations of Sutherland & Dopita (1993) and the cooling curves of 
Ferrara (1996, private communication) are adopted. For these high temperatures the cooling is dominated 
by thermal bremsstrahlung and by atomic hydrogen around 10 K. Unlike other treatments (Chiosi et 
al. 1997) we do not use cooling functions below 10 4 K, but solve the statistical equilibrium and the radiative 
transfer problem explicity (Spaans & van Dishoeck 1997; SN97). 

The fine-structure lines of C + , C, O, Si + , Si, S + , S, Fe + , and Fe are included, as well as the rotational 
lines of CO and H2. The vibrational excitation of the latter molecule is also included. The thermal balance 
of the molecular phase strongly depends on the chemistry and thereby influences the effective SF rate of 
the interstellar gas. A detailed treatment of optical depth effects in the cooling lines is therefore required 
to determine the various densities and temperatures of the components of the multi- phase ISM (NS97). 
Note that in rcdshift dependent calculations, the temperature of the Cosmic Microwave Background 
(CMB) should be taken explicitly into account for redshifts larger than 5. The latter effect is of particular 
importance for the excitation of CO molecules, and influences the effective cooling rate of molecular gas 
(Silk & Spaans 1997). Still, we do not expect the effective star formation rate to change much and the 
application of our results is likely to be to lower redshift objects. Therefore we have considered only a 3 K 
CMB. 



3. Results and Discussion 



3.1. Model Characteristics 
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In the computation of the star formation rates, the total body of gas simulated is M g = 10 10 M , 
i.e. the gas mass of the Milky Way. The effective star formation rates presented here should be scaled 
linearly when implemented in dynamical simulations with some mass resolution element M ros . To model 
the effect of outflows on the star formation process the body of gas is assumed to be embedded in a 
M pot = 5 x 10 11 M Q spherical potential with a r~ 2 density profile outside a constant density core (Spaans 
& Carollo 1997). The energy estimates of SN97 are used to determine when the gas becomes unbound and 
star formation is quenched. We have assumed a generic SN kinetic energy of 10 51 ergs, and that 10% of 
the initial SN energy ends up as turbulent kinetic energy of the interstellar clouds. The amount of kinetic 
energy transferred from the SN expansion front to the ISM is not that well determined. Nevertheless, recent 
one-dimensional numerical simulations by Thornton et al. (1997) indicate that 10% is a reasonable value. 
The value of M pot determines the binding energy of the gas in a linear fashion. The maximum possible star 
formation rate, i.e. when a galactic wind does not yet develop, therefore scales with M pot in our formalism. 
This scaling does not take into account the nature of the outflow which develops for a galaxy with a specific 
density profile. If the potential is not static the local density changes should be well captured by our 
formalism (barring shocks). The maximum star formation rate on the other hand would change, and our 
method does not accurately capture these effects. 

The characteristic time scale and absolute rate of star formation depend on the thermal and chemical 
balance of the ISM. For the numerical simulations presented here, a grid was set up in parameter space 
which includes total hydrogen number densities nn = 10~ 2 — 10 6 cm~ 3 and metallicities Z = 5 x 10~ 3 — 3 
Zq . The locally produced stellar radiation field and the input of kinetic and thermal energy by supernova 
explosions render our equation for the local star formation rate non-linear due to feedback effects. The 
numerical simulations have therefore been performed for a range of "background" SF rates Sb = 0.3 — 300 
M yr _1 . The SF rates at the low end of the background range can be much larger than 300 M Q yr _1 
initially, but as S increases, feedback effects in the form of outflows limit the sustained SF rate to ~ 10 3 
M yr _1 . This result depends on the adopted values of M g and M pot as discussed above. In practical 
applications, the background star formation rate should be defined be neighboring grid cells or a previous 
time step. 



3.2. Physical Processes 

Before turning to the specific results, it is good to bear in mind the following physical processes which 
determine to a large extent the conclusions of this paper. For metallicities of less than a percent of Solar, 
the chemistry is dominated by H, He, Li, their ions, H , and molecules like H 2 and HD. The latter two, 
together with atomic hydrogen, provide the cooling of the ambient gas. The dominant formation route of 
H2 involves the H~ radical rather than dust grains. A consequence, of importance for the results presented 
here, is that star formation is a global process. That is, the lack of dust allows stellar UV photons to 
traverse throughout the medium almost unattenuated and to influence the physical structure of the ISM 
over very large distances. In NS97 it is shown that this mode of star formation is self-regulating due to the 
strong dependence of the H~ abundance on the radiation field strength shortward of 912 A, and the fact 
that H 2 is only dissociated by photons longward of 912 A. The typical H 2 driven star formation rate for 
a proto-galactic disk is small and of the order of 0.01-0.1 M Q yr _1 . The relevance of this number for low 
surface brightness galaxies will be discussed in a forthcoming paper (Spaans, Mihos, & McCaugh 1998; in 
preparation) . 
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For metallicities in excess of a few percent of Solar, the above picture changes dramatically. The 
presence of metallic atoms enhances the cooling rate and allows the formation of molecules like OH, H2O 
and CO. Also the formation of H 2 proceeds much more efficiently on the surfaces of dust grains. Due to the 
increased dust opacity, the formation of stars is a local process. Stellar photons are absorbed close to the 
source and the presence of metallic atoms and molecules allows a large fraction of this input energy to be 
radiated away efficiently. 

The following phases can then be distinguished in the evolution of the ambient medium: Only a small 
fraction of baryonic matter is converted into stars due to H2 and HD cooling in a very low metallicity 
medium. The regions where the H 2 abundance is large are characterized by a kinetic temperature of ~ 1000 
K, due to vibrational line cooling, and an ambient pressure of 2 x 10 3 K cm~ 3 . The efficient formation 
of stars begins with a transition to a multi-phase ISM at 5 x 10 3 K cm~ 3 , which includes cold molecular 
clouds. This transition occurs at a metallicity of the order of ^0.03Z©. The multi-phase ISM is able to cool 
away most of the UV stellar photons through infrared lines and dust continuum emission. 



3.3. Model Results 

Figure 1 presents the general dependence on metallicity and hydrogen number density (for a mass 
M g ) and a background star formation rate 5b = 3 M Q yr _1 . The dots in Figure 1 (and Figure 2) indicate 
individual entries in Table 1. Note here that the model results presented in Table 1 have been interpolated 
with a low order polynomial to create a grid which is homogeneous in both density and metallicity. Around 
the threshold metallicity of ~ 0.03 Z©, the sustained SF rate is a strong function of the ambient enrichment. 
Because the cooling rate cannot keep up with the ambient heating, a small rise in the SF rate raises the 
temperature of the ISM to above ~ 1000 K where H 2 is collisionally destroyed. Conversely, a small increase 
in metallicity increases both the cooling and the H2 formation rate on dust grains, and therefore the 
sustained SF rate. Once a multi-phase ISM has been firmly established for metallicities roughly within an 
order of magnitude of Solar, the star formation rate does not vary with metallicity by more than a factor of 
two. This relative insensitivity to metallicity in a multi-phase ISM is due to the ability of the gas to cool 
through infrared line transitions and dust emission. 

Note that the 0.01 and 0.02Z© curves are the same for number densities below 0.03 cm -3 , but not 
the same as the 0.005 and 0.05Z© curves below that number density. Below a metallicity of ^O.OIZ©, a 
multi-phase ISM does not exist. The 0.005Z© curve therefore pertains to a medium cooled predominantly 
by H, H 2 and HD. Conversely, for metallicities close to the phase transition, but at low densities, the 
pressure-density curve which determines the phase structure of the ambient gas, does not vary much 
with metallicity. This is caused by the fact that for these parameters the heating and cooling rates vary 
proportionally. 

Figures 1 and 2 also illustrate the importance of feedback to limit the global SF rate to no more than 
<~ 10 3 M© yr _1 for the masses and ambient densities considered here. A system with (M pot = 5x 10 11 M©, 
M g = 10 10 M©) starts to unbound its ISM kinematically for a star formation rate of ~ 10 3 M© yr _1 , even 
though the thermal energy contributed to the molecular ISM can be largely radiated away through line and 
dust emission. This causes a break in the star formation rate with density for a given metallicity. For low 
densities and SF rates the slopes in Figures 1 and 2 are for a large part dictated by the assumed Schmidt 
law. This limit to the star formation rate should be indicative for the formation of L*-type galaxies, 
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and appears to be in agreement with current results for the Hubble Deep Field (Madau et al. 1996). A 
dynamical implementation of our results will yield a more accurate determination of these effects since they 
could account for re-accretion and the extent of the dark matter potential. 

It is good to take note of the time scales to reach chemical equilibrium in media of varying metallicity. 
For a low metallicity environment of 0.5 percent of Solar, the time to reach chemical equilibrium for the 
H2 abundance in a cloud of 100 cm -3 is approximately 500 Myr for a background radiation field of 0.1 in 
units of J_2i- In the multi-phase ISM the corresponding time scale is a bit less than 20 Myr. It is the 
H2 abundance which provides the initial cooling and drives the ion-molecule chemistry required for the 
formation of molecules like OH, CH, H 2 0, CO etc. The decrease in this time scale with metallicity therefore 
gives a direct measure of the efficiency with which stars can be formed. 

These time scales are also relevant to the blow-out phenomenon in flattened galaxies. A supernova 
shock front originating from the disk of a galaxy and travelling at 100 km s _1 will reach a height of 1 kpc 
in 10 Myr. This time scale is of the same order as the chemical equilibrium time as well as the cooling time 
of the gas. The hot metal-rich cavity gas therefore does not have enough time to cool and form (molecular) 
clouds which are more robust and could help to retain the metals, even in a multi-phase ISM (Spaans & 
Carollo 1997). 



3.4. Sensitivity to Model Assumptions and Parameters 

Some of the numbers which go into the calculation have a large influence on the final result and are 
also poorly determined. These will be discussed below. 

1) An optimal metallicity range of <~ 0.04 — 0.2Zq seems to exist, for which SF can proceed more 
efficiently at high densities (SN97) . This effect is due to the absence of magnetic support for low metallicity 
molecular clouds. A low ionization fraction yields an ambipolar diffusion time scale which is shorter than 
the free-fall time of a molecular cloud. The calculations presented here limit this increase in the SF rate 
to about of factor of two due to the grain chemistry adopted at high extinction. Deep inside molecular 
clouds grains are predominantly neutral and provide an important sine for the free electrons followed by 
rapid recombination with (metallic) ions (Lepp & Dalgarno 1988). This leads to a medium with a higher 
fractional ionization compared to the case where electron recombination dominates. The rates for these 
processes arc typically 10 -6 cm 3 s , but are not accurately known due to their dependence on the grain 
size distribution. Additional calculations at a selected number of densities which involve a chemical network 
with more detailed charge transfer effects and exchange rates larger by a factor of three, indicate that the 
star formation rates in the optimal metallicity range may be larger by about 100%. The more conservative 
estimates are retained here because only cloud collapse models could provide a definite answer. These seem 
to indicate that diffusion of the magnetic field aids the formation of stars, but is not a sufficient condition 
by itself (Gammie 1997, private communication). 

2) Related to the issue of molecular clouds at high visual extinctions, is the matter of grid resolution. 
The implemented method allows the gas to be resolved down to the 20 pc level, i.e. individual molecular 
clouds. Any increase in grid resolution is therefore not expected to change our results. Still, the star 
formation process itself involves the fragmentation of molecular cloud cores and smaller scales should 
therefore be probed to put the adopted Schmidt star formation law on a firmer basis. Such an increase in 
resolution must be accompanied by a hydrodynamic treatment of the gas to accurately model for instance 
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the effects of shocks. The latter are very dissipative phenomena and strongly influence the thermal balance 
of the ambient gas. It is expected that in a high metallicity medium, molecules will be (re-)formed efficiently 
in the shock's wake. Therefore, the energy will be efficiently radiated away in the form of long wavelength 
photons and the gas can condense to high densities. This will increase the effective star formation rate, 
although the magnitude of the effect lies beyond the scope of the present work. 

3) The heating rate due to supernova explosions can become comparable to the UV heating rate in 
some parts of the model galaxy. Nevertheless, the former number is not accurately known and we have 
therefore run models where the thermalized fraction of the input supernova energy is a factor of 5 larger. 
For low mctallicities in a H 2 cooling dominated medium, the temperature of the gas rises with a factor 
of a few, from ~ 1000 to ~ 4000 K. This increases the collisional destruction rate of molecular hydrogen, 
the dominant coolant, and therefore inhibits further star formation. Once a multi-phase ISM has been 
established, the ambient gas is capable of radiating away the extra input energy efficiently and the results 
for the star formation rates do not change by more than 30-50%. 

4) We have ignored the contribution from Type lb supernovae. It appears that radiative stellar winds 
from massive stars may be significantly suppressed at low metallicities, but that winds driven by binary 
star evolution could increase the number of Type lb supernovae in the early universe (Timmes 1998, private 
communication). Such an effect could increase the relative abundance of oxygen. Since atomic oxygen is 
an important coolant and strongly influences the chemistry of OH and CO, our results potentially can be 
altered by a higher Type lb supernova rate. For a higher oxygen production rate of a factor of three, it 
turns out that the cooling rate due to [OI] is not strongly influenced because the main cooling line, [OI] 63 
/tm, is optically thick. The abundance of CO is favorably increased and the molecular phase can cool more 
efficiently For the adopted factor of three, the star formation rate maximally increases by 40%. 

If one combines the uncertainties 1-4, it appears that the magnitude of the star formation rates 
presented here can change by at most a factor of four, and are typically accurate to a factor of two. We feel 
that the existence of a threshold metallicity, as depicted in Figures 1 and 2, is quite robust since it derives 
from a detailed treatment of the thermal balance of the ISM for a range of background radiation fields. 



4. Conclusions 

We have derived star formation rates in tabular form as functions of the ambient density, metallicity, 
and stellar feedback in proto-galaxies with a given total gas mass. We have explicitely included the effects 
of supernova explosions on the thermal balance of the ISM, and the transition to a multi-phase ISM. 
The strength of the method used to produce the effective star formation rates is the detailed treatment 
of molecular line cooling, chemistry, and thermal balance. The tabulated star formation rates can be 
implemented in N-body codes for the computation of star formation processes in merging galaxies and 
cosmological simulations. The accuracy of the presented numbers is expected to be of the order of a factor 
of two, and in any case the rates should not change by more than a factor of four. The full set of numerical 
results is included in Table 1, and can be obtained in machine readable form from MS. 

Although we have included the effects of feedback self-consistently, the present calculations cannot 
capture the effects of massive shocks associated with a large collection of SN explosions propagating through 
the ISM. A complete treatment of the latter requires the incorporation of the numerical scheme adopted 
here in a full hydro treatment. Investigations along these lines are being pursued. 
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Table 1 

Effective Star Formation Rates (SF)"ai^Functions of Metallicity, 
Background Star Formation Rate (BSF) and Density 
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a In M Q yr _1 . 

6 The notation aEb denotes a x 10 h . 
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Fig. 1. — The star formation rate as a function of hydrogen number density, for various metallicities, a total 
mass of M g and a background star formation rate of Sb — 3 M Q yr _1 . The various curves are labelled by 
the ambient metallicity relative to Solar. 



